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FINAL  TECHNICAL  REPORT 
FOR  AFOSR  GRANT  AFOSR  84-0113 
Principal  Investigators:  Henry  W.  Block,  Allan  R.  Sampson,  Thomas  H.  Savits 

1.  SUMMARY  OF  RESEARCH 

In  this  section  we  summarize  our  research  accomplishments  under  this  grant.  Our  research 
achievements  can  be  categorized  conveniently  into  five  areas:  reliability,  dependence  concepts,  applica¬ 
tions  of  dependence  concepts,  time  series,  and  other  research. 

A.  Reliability 

Under  this  grant,  various  topics  in  reliability  theory  have  been  studied.  These  include  multivariate 
nonparametric  classes,  maintenance  models  and  multistate  structures. 

In  technical  report  85-01,  recent  research  in  multivariate  nonparametric  life  classes  in  reliability  is 
surveyed.  The  multivariate  aging  notions  of  IFR,  IFRA,  NBU,  DMRL,  and  HNBUE  are  discussed  and 
some  of  their  relationships  and  properties  are  presented.  It  is  shown  in  technical  report  84-05  that  the 
IFRA  scaled-mins  class  is  closed  under  convolution.  A  multivariate  distribution  function  of  n  com¬ 
ponents  subjected  to  shocks  from  a  nonhomogeneous  Poisson  process  is  considered  in  report  85-07, 

It  is  shown  that  if  the  disu  ution  function  of  the  time  until  the  first  shock  is  IFR,  IFRA,  or  NBU,  then 
the  multivariate  lifelengths  of  the  n  components  belongs  to  a  corresponding  multivariate  IFR,  IFRA,  or 
NBU  class. 

A  survey  of  the  recent  research  concerning  maintenance  policies  is  presented  in  technical  report 
84-04.  Some  stochastic  comparisons  on  the  number  of  failures  and  removals  in  [0,s]  for  a  complete 
repair  model  under  age  and  block  replacements  are  obtained  in  86-08.  These  results  are  extended  for 
the  block  policy  case  in  technical  report  88-06  in  which  stochastic  comparisons  of  entire  processes  are 
obtained  for  the  complete  repair  model  as  well  as  for  the  minimal  repair  model.  The  technical  reports 
87-06  and  88-01  consider  a  very  general  cost  structure  for  a  maintained  system.  The  first  report 
obtains  a  cost  relationship  between  the  age  and  block  policies;  the  second  extends  the  relationship  to 
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include  discounting. 

In  technical  report  86-06  characterizations  of  nonhomogeneous  Poisson  processes  in  terms  of  cum- 
muiants  and  increments  of  the  counting  process  are  obtained.  An  estimate  of  the  stable  availability  for 
a  series  system  is  considered  in  87-02.  Some  related  asymptotic  normality  and  rates  of  convergence  are 
also  discussed.  Technical  report  87-03  investigates  a  class  of  multistate  systems  which  are  more 
series-like  or  more  parallel-like.  Some  closure  properties  are  presented  as  well  as  some  characteriza¬ 
tion  results.  Some  improved  reliability  bounds  for  systems  with  dependent  components  are  derived  in 
technical  report  87-10. 

B.  Dependence  Concepts 

Under  the  grant,  concepts  of  positive  dependence,  negative  dependence,  monotone  correlation  and 
some  applications  have  been  studied. 

Technica.  report  84-01  presents  a  unified  framework  for  studying  and  relating  positive  dependence 
orderings,  positive  dependence  properties  and  measures  of  positive  dependence.  The  latter  two  new 
concepts  are  formally  defined,  their  properties  are  discussed  and  interrelationships  among  these  three 
concepts  are  given.  In  technical  report  85-05  the  concepts  of  conditionally  more  positively  quadrant 
dependent  and  conditionally  more  dispersed  are  introduced  and  studied.  The  report  also  considers 
other  conditional  positive  dependence  concepts.  In  technical  report  85-lOa  representations  of  upper 
sets  contained  in  a  finite  two  or  three  dimensional  discrete  lattice  are  presented.  These  representations 
are  then  used  for  a  variety  of  applications  in  multivariate  positive  dependence.  In  technical  report  86 
-07,  the  nonparametric  maximum  likelihood  estimators  of  two  multivariate  c.d.f.’s  F  and  G  are 
obtained  under  the  constraint  that  F  is  stochastically  larger  than  G.  New  higher  order  two  and  three 
dimensional  product-type  lower  bounds  are  obtained  for  certain  multivariate  probabilities  in  technical 
report  88-08.  The  positive  dependence  conditions  sufficient  for  these  bounds  to  hold  are  described  and 
exemplified  by  the  multivariate  normal  and  Marshall-Olkin  exponential  distributions.  In  technical 
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report  87-05,  »  variety  if  <>evv  dependence  notions  arc  given  to  describe  the  dependence  relationships 
between  two  or  more  sets  of  variables.  Under  a  variety  of  different  notions  of  setwise  dependence, 
various  characterizations  of  setwise  independence  are  described  in  technical  report  87-04, 

Technical  report  84-02  studies  random  vectors  conditioned  on  events  involving  their  order  statis¬ 
tics.  The  conditioning  events  considered  occur  naturally  in  reliability  theory.  It  is  shown  that  these 
conditional  random  vectors  are  negatively  dependent  in  a  very  strong  sense.  In  technical  report  84-03 
a  concept  of  negative  dependence  by  stochastic  ordering  is  introduced.  It  is  shown  that  three  models 
of  negative  dependence  satisfy  it,  and  that  it  is  satisfied  by  some  frequently  used  multivariate  d.f.s. 
Technical  report  85-05  considers  conditional  negative  dependence  concepts  and  -<srr,*  nf  tb»,jr  onniira. 
tions. 

In  technical  report  85-10b  a  generalization  of  the  Hoeffding-formula  for  multivariate  cumulants  is 
presented.  This  result  is  then  used  to  characterize  independence.  Relationships  with  various  types  of 
dependence  are  also  given.  In  technical  reports  86-03  and  86-04,  a  bivariate  random  vector  assuming 
values  in  a  bivariate  lattice  is  considered.  The  probabilities  of  this  bivariate  random  vector  are 
presented  as  a  matrix  and  a  new  matrix  decomposition  of  Q  is  obtained.  This  decomposition  is  applied 
to  various  concepts  in  positive  dependence.  Furthermore,  a  number  of  results  concerning  the  concor¬ 
dant  monotone  correlation  and  the  sup-correlation  are  presented.  Additionally  considered  are  applica¬ 
tions  of  these  ideas  to  Hotelling’s  canonical  correlation  and  to  other  monotone  correlations. 


A  number  of  research  efforts  concerned  various  -,siuvs  dependence  notions  for  bivariate  and, 
more  generally,  multivariate  empirical  c.d.f.’s.  Through  invariance  arguments,  it  was  shown  that  these 
notions  can  be  studied  on  Sn,  the  set  of  all  permutations  of  { 1,  ...» n}.  In  technical  report  87-09,  a 

j 

variety  of  partial  orderings  on  Sn  are  considered  and  shown  to  be  related  to  certain  posith-e  depen-  - — 

dence  orderings  on  bivariate  c.d.f.’s.  The  concepts  are  further  amplified  upon  in  technical  report 


88-10.  In  technical  report  88-09,  certain  subsets  of  S„  =  Sn  x  Sn,  and  SnK  are  characterized  and  shown 


to  correspond  to  certain  multivariate  positive  dependence  concepts.  A  description  and  characterization 
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of  some  of  the  partial  orderings  on  Sn  by  means  of  metrics  on  Sn  is  given  in  technical  report  87-11. 

C.  Applications  of  Dependence  Concepts 

One  of  the  focuses  of  the  research  in  our  grant  concerning  dependence  is  on  the  applicability  of 
these  notions  to  statistical  and  data  analytical  issues. 

Technical  report  85-03  employs  positive  dependence  concepts  and  orderings  to  compare  and  evalu¬ 
ate  four  simultaneous  confidence  interval  procedures  for  estimable  functions  in  the  general  linear 
model.  The  probability  inequalities  of  technical  report  88-08  are  used  in  technical  report  88-07  to 
develop  an  "optimized"  product-type  probability  bound  which  is  important  in  a  variety  of  applications, 
e.g.,  change  point  detection  and  group  sequential  trials.  Under  commonly  occuring  positive  depen¬ 
dence  structures,  this  new  bound  is  shown  to  be  superior  to  previously  obtained  optimized  Bonferroni- 
type  bounds  of  Hunter  and  Worsley.  In  technical  report  87-08,  some  recent  multivariate  probability 
inequalities  are  used  to  obtain  simultaneous  inference  procedures  consistent  with  the  Benard-van 
Elteren  nonparametric  procedure,  which  is  used  to  compare  multiple  treatments,  in  the  presence  of 
blocking  effects  and  possibly  unequal  sample  sizes.  Technical  report  85-09  utilizes  notions  of  mono¬ 
tone  correlation  to  develop  a  methodology  to  merge  in  a  unidimensional  scale,  judgements  of  two 
evaluators  whose  opinions  (possibly  on  dissimilar  scales)  about  N  objects  are  given  in  bivariate  ordinal 
contingency  table  form. 

Two  other  research  efforts  focus  on  the  application  of  the  notions  of  association  and  positive  qua¬ 
drant  dependence  to  ordinal  contingency  tables  and  bivariate  empirical  distributions.  In  technical  report 
88-03,  an  efficient  computational  algorithm  is  given  for  checking  association  in  bivariate  ordinal  con¬ 
tingency  tables.  For  even  moderate  sized  tables,  this  algorithm  is  shown  to  reduce  the  computational 
effort  by  potentially  several  orders  of  magnitude.  Technical  report  88-04  gives  a  new  characterization 
for  bivariate  empirical  distributions  to  be  positively  quadrant  dependent.  Based  upon  this  characteriza¬ 
tion,  efficient  graph-theoretic  algorithms  are  given  to  enumerate  and  generate  all  positively  quadrant 
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dependent  empirical  distributions  of  a  fixed  sample  size. 

D.  Time  Series 

Another  area  of  research  for  this  grant  has  been  times  series  with  particular  attention  to  modeling 
nonnormal  multivariate  processes. 

These  accomplishments  in  time  series  are  given  in  four  technical  reports.  In  technical  reports 
85-02  and  86-01  various  time  series  with  bivariate  exponential  and  geometric  marginals  are  presented. 

It  is  shown  that  dependency  on  the  past  is  relatively  simple  and  explicit.  The  theory  of  positive  depen¬ 
dence  is  used  to  obtain  association  between  random  vectors.  Special  cases  of  these  processes  are  sta¬ 
tionary  with  well  known  bivariate  exponential  and  geometric  distributions.  A  similar  project  airing 
Dirichlet  distributions  is  reported  on  in  report  87-01,  In  87-07  a  more  general  case  with  arbitrary  con¬ 
tinuous  marginals  is  studied. 

E.  Other  Research 

A  number  of  other  various  research  topics  were  investigated,  in  part,  by  co-researchers  and  visiting 
researchers  who  participated  in  our  grant.  These  are  now  summarized. 

In  technical  report  85-04  log-linear  models  are  applied  to  ordinal  longitudinal  data,  which  were 
collected  to  study  the  mental  health  effects  of  the  Three  Mile  Island  nuclear  accident.  Technical  report 
85-06  proposes  a  new  method  for  assessing  degree  of  agreement  among  raters,  and  discusses  the 
advantages  of  this  method.  In  technical  report  85-08  a  methodology  for  determining  standards  based 
on  multivariate  tolerance  regions  is  presented.  Their  use  in  quality  control  is  discussed.  In  technical 
report  85-1 1  the  question  of  the  effect  of  adenoidectomy  on  the  bearing  loss  of  children  with  SOM  is 
addressed  from  the  statistical  point  of  view,  and  in  technical  report  85-12  the  issue  of  evaluation  of 
scientific  information  is  addressed  using  statistical  methodology.  Technical  report  85-13  presents  the 
cross  prediction  problem  in  a  general  linear  regression  model.  Generalizations  of  this  problem  are  dis¬ 
cussed  and  applications  are  illustrated.  In  technical  report  85-14  two  characterizations  for  linear 
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exponential  families  are  presented.  Technical  report  86-02  analyzes  randomized  cloud  seeding  data. 
Technical  report  86-05  presents  an  alternative  method  for  obtaining  M.L.E  estimates  in  the  discrete 
time  mover-stayer  model,  and  illustrates  the  resuits  with  a  data  applicanon.  A  review  of  Multiple 
Interpretation  of  Clinical  Laboratory  Data  by  A.  Albert  and  E.  K.  Harris  (1987)  appears  in  report 
87-12.  Technical  report  87-02  gives  exact  power  and  sample  size  calculations  for  planning  and  design 
of  observational  studies.  In  technical  report  88-05  the  reliability  of  a  cytolytic  assay  performed  at  the 
Pittsburgh  Cancer  Institute  is  described  and  discussed. 

2.  VISITORS 

During  the  grant  period  May  1984  to  present,  a  large  nvmber  of  researchers  have  visited  for  vari¬ 
ous  time  periods.  A  list  of  these  visiting  scholars  with  their  affiliation  and  visit  dates  is  as  follows: 


M.  Abdel-Hammed  (University  of  Kuwait) 

M.  Bhattacharjee  (Indian  Institute  Management) 
P.  Boland  (University  College,  Dublin) 

J.  Bryant  (University  of  Cincinnati) 

S.  Cambanis  (University  of  North  Carolina) 

D.  Chhetry  (Katmandu  University) 

M.  L.  Eaton  (University  of  Minnesota) 

E.  El-Neweihi  (University  of  Illinois,  Chicago) 
E.  Enns  (University  of  Calgary) 

Z.  Fang  (University  of  Pittsburgh) 

C.  Fuchs  (Tel-Aviv  University) 

W.  Griffith  (University  of  Kentucky) 


August  2-11,  1984 

July  2-4,  1984 

June  9-15,  1985 

September  1987  -  August  1988 

March  31,  1988 

September  1,  1985-Augst  31,  1987 

February  5,  1988 

July  2-4,  1984;  May  13-17,  1986 

October  10-11,  1985 

September  1,  1986-August  31,  1987 

September  1,  1985-May  31,  1986 

July  22- August  1,  1984; 

July  10-12,  15-17,  1985; 

May  19-22,  1986 


G.  Kimeldorf  (University  of  Texas,  Dallas) 


July  18-20,  1984 
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N.  Langberg  (Haifa  University)  June  15- August  31,  1984; 

February  1 -March  1,  1985; 

June  1-July  1,  1985; 

February  1-March  1,  1986; 
September  1,  1986- August  31,  1987 
May  -  July  1988 

M.-L.  Lee  (Boston  University)  June  24-July  7,  1984; 

July  21-26,  1985 

S.  May  (Camegie-Mellon  University)  September  1,  1985-April  30,  1986 

T.  Nguyen  (Bowling  Green  University)  June  26-July  7,  1986 

B.  Schweitzer  (University  of  Massachusetts)  October  1987 

M.  Shaked  (University  of  Arizona)  December  7,  1987 

N.  Singpurwalla  (George  Washington  University)  November  14-16,  1985 

E.  Sungar  (Carnegie-Mellon  University)  January  1-May  31,  1986 

J.  Tiago  de  Oliveira  (Academy  of  Sciences  of  Lisbon)  July  1,  1987 

G.  Tsaturyan  (Florida  State  University)  April  1988 

L.  Whitaker  (University  of  CA,  Santa  Barbara)  July  1,  1985-March  31,  1986 

July  11-16,  1986 


3.  TITLES  AND  ABSTRACTS  OF  TECHNICAL  REPORTS 


REPORT  NO. 


TITLE 


AUTHOR(S) 


84-01  A  FRAMEWORK  FOR  POSITIVE  DEPENDENCE  George  Kimeldorf 

Allan  R.  Sampson 


Abstract  -  This  paper  presents  for  bivariate  distributions  a  unified  framework  for  studying  and  relat¬ 
ing  three  basic  concepts  of  positive  dependence.  These  three  concepts  are  positive  dependence  order¬ 
ings,  positive  dependence  properties  and  measures  of  positive  dependence.  The  latter  two  concepts  are 
formally  defined  and  their  properties  discussed.  Interralationships  among  these  three  concepts  are 
given,  and  numerous  examples  are  presented. 


84-02 


PROBABILITY  INEQUALITIES  VIA  NEGATIVE 
DEPENDENCE  FOR  RANDOM  VARIABLES 
CONDITIONED  ON  ORDER  STATISTICS 


Henry  W.  Block 
Vanderlei  Bueno 
Thomas  H.  Savits 
Moshe  Shaked 


Abstract  -  Distributions  are  studied  which  arise  by  considering  independent  and  identically  distributed 
random  variables  conditioned  on  events  involving  order  statistics.  It  is  shown  that  these  distributions 
are  nega-  tively  dependent  in  a  very  strong  sense.  Furthermore,  bounds  are  found  on  the  distribution 
functions.  The  conditioning  events  considered  occur  naturally  in  reliability  theory  as  the  time  to  sys¬ 
tem  failure  for  k-out-of-n  systems.  An  application  to  systems  found  with  "second-hand”  components  is 
given. 


84-03  A  CONCEPT  OF  NEGATIVE  DEPENDENCE  Henry  W.  Block 

USING  STOCHASTIC  ORDERING  Thomas  H.  Savits 

Moshe  Shaked 


Abstract  -  A  concept  of  negative  dependence  called  negative  dependence  by  stochastic  ordering  is 
introduced.  This  concept  satisfies  various  closure  properties.  It  is  shown  that  three  models  for  nega¬ 
tive  dependence  satisfy  it  and  that  it  implies  the  basic  negative  orthant  inequalities.  This  concept  is 
also  satisfied  by  the  multinomial,  multivariate  hypergeomeu»c,  Diiuduet  and  DuiddU  compound  multi¬ 
nomial  distributions.  Furthermore,  the  joint  distribution  of  ranks  of  a  sample  and  the  multivariate  nor¬ 
mal  with  non-  positive  pairwise  correlations  also  satisfy  this  condition.  The  positive  dependence  ana¬ 
log  of  this  condition  is  also  studied. 
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84-04  PREVENTIVE  MAINTENANCE  POLICIES  Henry  YV.  Block 

YVagner  S.  Borges 
Thomas  H.  Savits 

Abstract  -  Recent  research  concerning  maintenance  policies  with  time  dependent  cost  and  probabili¬ 
ties  is  surveyed. 


84-05  CONVOLUTION  OF  THE  IFRA  SCALED-MINS  CLASS  Emad  El-Neweihi 

Thomas  H.  Savits 


Abstract  -  The  class  of  nonnegative  random  vectors  T  =  (Tj,  .  .  .  ,Tn)  for  which  min  a,Tj  is 

1  l  i  s  n 

IFRA  for  all  0  <  2lx  <  «>,  i  =  1,  .  .  .  ,n,  is  closed  under  convolution. 


85-01 


MULTIVARIATE  NONPARAMETRIC  CLASSES 
IN  RELIABILITY 


Henry  W.  Block 
Thomas  H.  Savits 


Abstract  -  Recent  research  in  multivariate  nonparametric  classes  in  reliability  is  surveyed. 
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85-02  MOVING  AVERAGE  MODELS  WITH  BIVARIATE  Naftali  A.  Langberg 

EXPONENTIAL  AND  GEOMETRIC  DISTRIBUTION  David  S.  Stoffer 

AbsL  act  -  Two  classes  of  finite  and  infinite  moving  average  sequences  of  bivariate  random  vectors 
are  considered.  The  first  class  has  bivariate  exponential  marginals  while  the  second  class  has  bivariate 
geometric  marginals.  The  theory  of  positive  dependence  is  used  to  show  that  in  vanous  cases  the  two 
classes  consist  of  associated  random  variables.  Association  is  then  applied  to  establish  moment  ine¬ 
qualities  and  to  obtain  approximations  to  some  joint  probabilities  of  the  bivariate  processes. 


85-03  SIMULTANEOUS  CONFIDENCE  INTERVALS  Camil  Fuchs 

FOR  THE  GENERAL  LINEAR  MODEL  Allan  R.  Sampson 

Abstract  -  This  paper  describes,  compares  and  illustrates  four  known  simultaneous  confidence  interval 
procedures  for  the  general  linear  model.  It  is  shown  that  three  of  the  four  procedures  can  be  ordered 
by  the  increasing  length  of  their  confidence  intervals  as  follows:  (I)  Sidak  multivariate  t-distnbution 
intervals,  (II)  Sidak  independent  t-distribution  intervals,  and  (III)  Bonferroni  intervals.  Tnere  is  no  uni¬ 
form  ordering  between  the  confidence  intervals  yielded  by  the  fourth  procedure  (SchefffeJ  and  the 
other  three;  these  relationships  are  explored  numerically. 


85-04  ORDINAL  LOG-LINEAR  MODELS  AND  THE  TMI  DATA  Sherryl  May 

Abstract  -  Log-linear  models  appropriate  for  data  arranged  in  ordinal  contingency  tables  have  been 
proposed  and  discussed  by  authors  such  as  Fienberg  (1982)  and  Agresti  (1983).  Such  models  are 
applied  to  ordinal  longitudinal  data  collected  on  the  subject  of  possible  mental  health  effects  of  the 
Three  Mile  Island  nuclear  accident.  The  results  are  discussed. 


85-05  CONDITIONALLY  ORDERED  DISTRIBUTIONS  H.nry  R.  Block 

Allan  R.  Sampson 


Abstract  -  The  concepts  of  conditionally  more  positively  quadrant  dependent,  and  conditionally  more 
dispersed  are  introduced  and  studied.  Based  on  these  two  concepts,  new  conditions  are  given  for  mul¬ 
tivariate  c.d.f.’s  F  and  G  so  that  Eph(X)  >  Ech(X)  for  suitable  h(X).  Special  cases  include  the  mul¬ 
tivariate  normal  distribution  and  elliptically  contoured  distributions.  Conditional  positive  and  negative 
dependence  concepts  as  well  as  applications  to  the  Farlie  Gumbel-Morgenstern  distribution  are  also 
considered. 
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85-06  MODELLING  INTER-RATER  AGREEMENT  Sherryl  May 

WITH  ORDINAL  LOG-LINEAR  MODELS 


Abstract  -  A  new  method  for  assessing  the  degree  of  agreement  among  raters  is  proposed.  This 
method  makes  use  of  the  recently-developed  methodology  of  ordinal  log-linear  models.  The  advan¬ 
tages  of  this  new  method  are  discussed  and  explained.  Several  examples  of  the  implementation  of  the 
method  are  then  presented. 


85-07  SOME  MULTIVARIATE  DISTRIBUTIONS  DERIVED  Thomas  H.  Savits 

FROM  A  NONFATAL  SHOCK  MODEL 


Abstract  -  A  nonhomogeneous  Poisson  shock  model  has  a  continuous  mean  function  A(t).  The  kth 
shock  Sk  causes  simultaneous  failuie  of  die  components  j  eJc{l,...,n}  with  probability  Pj(Sk). 

If  Tj  is  the  lifetime  of  component  j,  it  is  shown  that  (Tj,  .  .  .  >Tn)  belongs  to  various  multivariate  non- 
parametric  life  classes  depending  on  the  life  class  of  F(t)  =  exp{-A(t)}. 


85-08  MULTIVARIATE  TOLERANCE  REGIONS  AND  F-TEST  Camil  Fuchs 

R.  Kenett 


Abstract  -  The  determination  of  standards  or  norms  is  a  central  issue  in  quality  control  procedures. 

In  this  article  we  present  a  methodology  for  determining  standards  based  on  multivariate  tolerance 
regions.  We  study  the  correspondance  between  the  F-test  at  significance  level  a  and  tolerance  regions 
with  parameters  (r,5). 

The  methodology  and  statistical  issues  arising  from  its  application  are  illustrated  with  two  real  data 
problems.  In  both  cases  the  attributes  of  the  tested  elements  are  compared  with  tolerance  regions 
obtained  from  a  reference  sample.  In  the  first  example,  both  the  tested  elements  and  the  reference 
sample  originate  from  the  same  target  ("in  control")  population.  The  significance  levels  of  the  result¬ 
ing  tests  are  estimated  by  a  crossvalidation  procedure.  The  tested  elements  in  the  second  sample  are 
ffom  an  "out  of  control"  process.  The  percent  of  elements  detected  as  not  belonging  to  the  target 
population  constitutes,  in  that  case,  the  empirical  power  of  the  test. 
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85-09  MERGING  DISSIMILAR  ORDINAL  Jerrold  H.  May 

CATEGORICAL  JUDGMENTS  Allan  R.  Sampson 

Abstract  -  Two  judges  are  presented  with  N  objects  and  are  asked  separately  to  evaluate  each  object 
according  to  a  fairly  specific  qualitative  criterion.  Each  judge  is  allowed  to  develop  his  own  unique 
ordinal  categorical  scale.  Thus  the  data  consists  of  N  pairs  of  ordinal  categorical  data  on  dissimilar 
scales.  We  give  an  approach  which  merges  the  separate  bivariate  ordinal  scales  into  a  single  numerical 
value  providing  a  rank  order  of  all  N  objects.  To  do  so,  a  model  is  first  postulated,  and  its  parameters 
are  estimated  using  a  technique  involving  the  non-linear  minimization  of  a  certain  natural  squared  error 
optimality  criterion.  The  computational  aspects  of  this  optimization  are  discussed.  The  procedure  is 
applied  to  a  large  data  set  involving  financial  evaluation. 


85-10a  POSITIVE  DEPENDENCE.  UPPER  SETS,  Allan  R.  Sampson 

AND  MULTIDIMENSIONAL  PARTITIONS  Lyn  R.  Whitaker 

Abstract  -  Let  U  be  an  upper  set  contained  in  the  finite  discrete  lattice 

L  =  {1,  .  .  .  ,D(1)}  x  •  •  ■  x  {1,  .  .  .  ,D(p)}.  Representations  for  U  are  obtained  and  shown  to 
correspond  to  certain  multidimensional  partitions  of  integers.  It  is  shown  that  for  p=3,  the  number  of 

D(3) -l  fryn  )-HTX2)+i  1  flX  1  )+i  ] 

possible  upper  sets  in  L  is  J"[  ^  v  £)(1)  ^  \  [  D(l)  J‘  Various  other  representation  and 

enumeration  results  are  obtained  for  related  settings.  Considered  are  a  variety  of  applications  to  mul¬ 
tivariate  positive  dependence  and  various  notions  in  statistics  and  reliability  theory. 


85-10b 


A  MULTIVARIATE  EXTENSION 
OF  HOEFFDING’S  LEMMA 


Henry  W.  Block 
Zhaoben  Fang 


Abstract  -  Hoeffding’s  Lemma  gives  an  integral  representation  of  the  covariance  of  two  random  vari¬ 
ables  in  terms  of  difference  between  their  joint  and  marginal  probability  functions,  i.e., 


Cov(X,Y)  =  JT  £  {P(X  >  x,  Y  >  y)  -  P(X  >  x)P(Y  >  y)}dxdy 


This  identity  has  been  found  to  be  a  useful  tool  in  studying  the  dependence  structure  of  various  ran¬ 
dom  vectors.  A  generalization  of  this  result  for  more  than  2  random  variables  is  given.  This 
involves  the  following  integral  representation  of  the  multivariate  joint  cummulart  for  r  >  1: 


Cum(Xl5 .  .  .  ,Xr)  =  J-J CumOcx,  (xi),  •  •  •  .Xxr(xr))dxr"dxr 


where  Xx,(xi)  =  1  if  X;  >  xi;  0  otherwise.  Applications  of  this  result  include  characterizations  of 
independence.  Relationships  with  various  type?  of  dependence  are  also  given. 
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85-11  ADENOIDECTOMY  VERSUS  VT  Camil  Fuchs 

Abstract  -  This  article  focuses  on  the  question  of  the  effect  of  adenoidectomy  on  the  hearing  loss  of 
children  with  SOM.  The  empirical  results  originate  from  a  recent  study  which  cc  nsiders  various  facets 
related  to  theory  that  correlates  adenoidectomy  with  the  cure  of  SOM  and  AOM.  The  question  of 
hearing  loss  is  one  of  the  five  hypotheses  addressed  in  that  study  Here  we  focus  on  the  models,  on 
the  assumptions  and  on  the  empirical  results  from  several  studies.  The  results  from  three  original  stu¬ 
dies  are  cross  validated  with  those  from  other  published  investigations. 


85-12  EV  ALUATION  OF  SCIENTIFIC  INFORMATION  Camil  Fuchs 

Abstract  -  This  paper  deals  with  issues  of  evaluation  of  scientific  information  from  a  statistical 
viewpoint  and  was  presented  at  the  International  Symposium  on  Acute  and  Secretory  Otitis  Media. 


85-13  AN  INTERPRETATIVE  FRAMEWORK  FOR  Allan  R.  Sampson 

REGRESSION  AND  LINEAR  PREDICTION  David  S.  Stoffer 

Abstract  -  When  (X,Y)'  N((|J.j,JJ.2)^  £),  X  positive  definite,  the  values  of  a,(3  which  minimize  the 
weighted  risk  E{w[Y-(a+[3X)]2  +  (l-w)[X-(Y  -a)/J3]2}  zre  obtained.  This  is  termed  the  cross¬ 
prediction  problem.  Generalizations  of  this  problem  are  presented  and  related  to  other  regression  and 
multivariate  techniques  such  as  principal  components,  thereby  providing  a  unified  framework  for  many 
of  these  concepts.  The  various  generalizations  are  examined  in  detail  and  illustrated  in  applications  to 
a  well  known  data  set 


85-14  ON  THE  LINEAR  EXPONENTIAL  FAMILY  Naftali  A.  Lang^ert 

WITH  A  QUADRATIC  VARIANCE  FUNCTION  Osnat  Strainer 

Abstract  -  It  is  shown  that  for  every  triple  (A,BJD),  A  >  0  and  B2  -  4AD  <  0  there  is  a  linear 
exponential  family  (LEF)  of  distributions  s.t.  its  variance,  as  a  function  of  the  mean  Ql),  is  given  by 
Aji2  +  B|l  +  D.  It’s  also  shown  that  if  for  a  triple  (A,B  JD),  A  <  0  and  B2  -  4AD  >  0  there  is  a 
LEF  with  variance  function  equal  to  Afj.2  +  Bfi  +  D  then  -A-1  is  a  positive  integer  and  the  LEF  is  a 
linear  transformation  of  a  binomial  LEF. 
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86-01  BIVARIATE  EXPONENTIAL  AND  GEOMETRIC 

AUTOREGRESSIVE  AND  AUTOREGRESSIVE 
MOVING  AVERAGE  MODELS 

Abstract  -  We  present  autoregressive  (AR)  and  autoregressive  moving  average  (ARMA)  processes 
with  bivariate  exponential  (BE)  and  bivariate  geometric  (BG)  distributions.  The  theory  of  positive 
dependence  is  used  to  show  that  in  various  cases,  the  BEAR,  BOAR,  BEARMA,  and  BGARMA 
models  consist  of  associated  random  variables.  We  discuss  special  cases  of  the  BEAR  and  BGAR 
processes  in  which  the  bivariate  processes  are  stationary  and  have  well  known  bivariate  exponential 
and  geometric  distributions. 


86-02  HISTORICAL  ANALYSIS  OF  THE  EFFECT  Yoav  Benjamini 

OF  CLOUD  SEEDING  ON  RUNOFFF:  Camil  Fuchs 

CONDITIONAL  VS  UNCONDITIONAL  ANALYSIS 

Abstract  -  Randomized  cloud  seeding  experiments  had  been  conducted  in  Israel  since  the  early  six¬ 
ties.  The  decision,  based  on  the  analysis  of  rainful  data,  was  to  seed  operationally  every  possible  day. 

In  this  article  we  present  a  method  to  study  the  effect  of  cloudseeding  on  the  runoff,  by  a  double 
regression  approach  that  uses  historical  data  on  both  flow  measurements  as  well  as  rainfall  from  target 
and  control  areas.  It  was  found  that  the  magnitude  of  the  standard  errors  of  the  parameter  which  meas¬ 
ures  the  effect  of  seeding  on  runoff,  depend  largely  on  whether  the  problem  is  viewed  conditionally  or 
unconditionally.  As  a  result,  the  seeding  effect  is  labeled  highly  significant  or  nonsignificant  depend¬ 
ing  on  the  conditionality  assumption.  Jackknife  and  bootstrap  estimators  of  the  standard  deviation 
were  called  for  aid  and  compared  with  crude  straightforward  estimators.  The  performance  of  the  vari¬ 
ous  estimates  of  the  standard  deviations,  in  our  situation,  was  studied  using  a  small  simulation  study. 
The  results  illustrate  the  contention  that  in  cases  in  which  the  investigator  is  uncertain  of  the  effects  of 
conditionality,  the  jackknife  and  the  bootstrap  present  a  simple  and  appealing  alternative. 


Henry  W.  Block 
Naftali  A.  Langberg 
David  S.  Stoffer 


86-03  A  PROJECTION  DECOMPOSITION  FOR  BIVARIATE  Allan  R.  Sampson 

DISCRETE  PROBABILITY  DISTRIBUTIONS  Devendra  Chhetry 

Abstract  -  Let  Q  =  (Prob(X  =  X;,Y  =  yj)}  for  Xj<  •  •  •  <xm,  yj<  •  •  •  <yn.  A  new  matrix  decom¬ 
position  of  Q  is  given  in  terms  of  certain  projections  on  linear  spaces  related  to  the  marginal  probabili¬ 
ties.  It  is  shown  that  this  decomposition  implies  Fisher’s  canonical  decomposition  and  also  a  represen¬ 
tation  important  in  positive  dependence.  Additionally  considered  are  applications  of  these  ideas  to  the 
concordant  monotone  correlation,  the  maximal  correlation,  and  Hotelling’s  canonical  correlation. 
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86-04 


MONOTONE  CORRELATION  AND 
MONOTONE  DISJUNCT  PIECES 


Devendra  Chhetry 
Jan  deLeeuw 
Allan  R.  Sampson 


Abstract  -  Suppose  X,Y  are  random  variables  taking  values  on  the  lattice 
{x!  <...<  xm}  x  {y j  <...<  y0}  with  Q  =  {Pro  )(X=x„  Y=yj)}.  Let  Pcmc(Q)  and  Pdmc(Q)  be  ^ 
concordant  and  discordant  monotone  correlations  defined,  respectively,  by  the  maximum  and  minimum 
of  (p(f(X),g(Y))}  over  all  f,g  increasing  with  nonzero  variances.  A  number  of  results  concerning 
Pcmc(Q)  and  Pdmc(Q)  and  their  evaluations  are  obtained.  One  result  shows  that  Pcmc(Q)  =  1.  if 
and  only  if  Q  consists  of  at  least  two  increasing  disjunct  pieces,  i.e.,  Q  =  Diag(Q],Q2).  Necessary 
and  sufficient  conditions  are  also  given  for  Pcmc(Q)  =  Pdmc(Q)- 


86-05  MAXIMUM  LIKELIHOOD  ESTIMATIONS  IN  Camil  Fuchs 

EXTENDED  MOVER-STAYER  MODELS 


Abstract  -  The  discrete  time  mover-stayer  model  is  a  useful  model  for  studying  change  over  time  in 
heterogeneous  populations.  In  this  report,  an  alternative  method  for  obtaining  maximum  likelihood 
estimates  of  the  parameters  of  the  mover-stayer  model  based  on  the  EM  algorithm  is  presented  and 
generalizations  of  the  basic  model  are  considered.  The  method  and  the  model  are  illustrated  with  data 
from  a  community  based  study  and  changes  in  mental  health  status  over  a  one  year  period. 


86-06  CHARACTERIZATION  OF  NONHOMOGENEOUS  Zhaoben  Fang 

POISSON  PROCESSES  VIA  MOMENT  CONDITIONS 


Abstract  -  The  property  of  independent  increments  is  one  of  the  most  important  for  defining  both  the 
homogeneous  and  nonhomogeneous  Poisson  process. 

In  this  paper  we  give  two  ways  to  relax  this  requirement  and  characterize  the  nonhomogeneous 
Poisson  process  by  some  moment  conditions. 

One  result  is  that  a  counting  process  {N(t),t  >  0}  with  finite  moments  of  all  orders  is  a  nonhomo¬ 
geneous  Poisson  process  with  mean  functions  m(t)  =  EN(t)  if  and  only  if  for  any  tj,  i  =  1,  .  .  .  ,k 

cum  (N(tt),  .  .  .  ,N(tk))  =  min  EN(tj) 

l<i<k 


where  cum  (•)  is  the  joint  multivariate  cumulant. 

A  second  result  is  that  if  increments  on  any  interval  are  Poisson  distributed  and  an  exchangeable 
condition  is  assumed  then  the  process  is  nonhomogeneous  Poisson.  This  extends  Renyi’s  (1967) 
result. 


86-07 


ESTIMATION  OF  MULTIVARIATE  DISTRIBUTIONS 
UNDER  STOCHASTIC  ORDERING 


Allan  R.  Sampson 
Lyn  R.  Whitaker 


Abstract  -  Let  F  and  G  be  the  cdfs  of  two  p-dimensional  multivariate  distribution,  such  that  F  is  sto¬ 
chastically  larger  than  G.  A  straightforward  derivation  is  given  of  a  the  generalized  maximum  likeli¬ 
hood  estimators  of  F  and  G  based  on  the  random  samples  from  each  population.  An  algorithmic 
approach  to  computing  these  estimators  is  described  and  motivating  numerical  examples  are  discussed. 
The  special  case  when  F  and  G  correspond  to  multivariate  ordinal  contingency  tables  is  also  presented. 
The  relationship  of  these  to  those  of  Robertson  and  Wright  (Ann.  Statist  2  (1974)  528-534)  is  con¬ 
sidered. 


86-08  COMPARISONS  OF  REPLACEMENT  POLICIES  Naftali  A.  Langberg 


Abstract  -  For  independent  random  lifelengths  of  the  units  in  use  stochastic  comparisons  of  the 
number  of  failures  and  removal  in  [0,s]  under  age  and  block  replacement  policies  are  performed.  A 
new  concept  of  NBUfNWU)  in  sequence  is  introduced. 


87-01 


A  NONLINEAR  AUTOREGRESSIVE  MODEL  IN  DIRICHLET 
RANDOM  VECTORS  -  THEORY,  ESTIMATION,  AND  PREDICTION 


N.  A.  Langberg 
D.  S.  Stoffer 


Abstract  -  We  construct  a  nonlinear  autoregressive-type  model  which  can  be  used  for  modeling  and 
forecasting  vector  processes  in  which  the  distribution  of  the  random  vector  at  each  point  in  time  or 
space  is  a  Dirichlet  distribution.  That  is,  the  model  is  useful  for  modeling  and  forecasting  vector 
processes  in  which  the  data  are  proportions  and  constrained  so  that  the  sum  of  the  proportions  are 
equal  to  one  at  each  point  in  time  or  space.  Estimators  of  the  parameters  of  the  model  are  obtained  via 
two  different  methods.  First  we  obtain  maximum  likelihood  estimators,  and  then  we  obtain  uniformly 
minimum  variance  unbiased  estimators  of  the  model  parameters.  Finally,  we  present  an  example  of  the 
kind  of  data  that  can  be  handled  by  the  model,  and  discuss  the  modeling,  estimation  and  prediction 
methods  for  the  given  data  set 
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87-02  AN  ESTIMATION  OF  AVAILABILITY  OF  SERIES  SYSTEM  Jie  Mi 

Abstract  -  Suppose  a  system  S  consists  of  n  s-independent  subsystems  Sj  (j  =  1 ,2,  .  .  .  ,  n)  con¬ 
nected  in  series.  Denote  the  lifetime  and  repair  time  of  subsystem  Sj  by  Uj  and  Dj,  respectively. 

Based  on  the  sample  of  complete  periods  for  each  subsystem  we  give  an  approximate  interval  estimate 
of  the  stable  availability  of  the  system  S  in  three  cases:  (i)  Each  Uj  and  Dj  has  exponential  distribu¬ 
tion;  (ii)  Each  Uj  is  distributed  exponentially  but  each  Dj  has  log-normal  distribution;  (iii)  Both  Uj  and 
Dj  are  distributed  according  to  log-normal  law.  Also  the  related  asymptotic  normality  and  the  rate  of 
convergence  are  discussed. 


87-03 


L-SUPERADDITIVE  STRUCTURE  FUNCTIONS 


Henry  W.  Block 
William  S.  Griffith 
Thomas  H.  Savits 


Abstract  -  Structure  functions  relate  the  level  of  operations  of  a  system  as  a  function  of  the  level  of 
the  operation  of  its  components.  In  this  paper  structure  functions  are  studied  which  have  an  intuitive 
property,  called  L-superadditive  (L-subadditive).  Such  functions  describe  whether  a  system  is  more 
series-like  or  more  parallel-like.  L-superadditive  functions  are  also  known  under  the  names  supermo- 
dular,  quasi-monotone  and  superadditive  and  have  been  studied  by  many  authors.  Basic  properties  of 
both  discrete  and  continuous  (i.e.,  taking  a  continuum  of  values)  L-superadditive  structure  functions  are 
studied.  For  binary  structure  functions  of  binary  values,  El-Neweihi  (1980)  showed  that  L- 
superadditive  structure  functions  must  be  series.  This  continues  to  hold  for  binary  valued  structure 
functions  even  if  the  component  values  are  continuous  (see  Proposition  3.1).  In  the  case  of  non-bina:y 
valued  structure  functions  this  is  no  longer  the  case.  We  consider  structure  functions  taking  discrete 
values  and  obtain  results  in  various  cases.  A  conjecture  concerning  the  general  case  is  made 


87-04 


SETWISE  INDEPENDENCE  FOR  SOME 
DEPENDENCE  STRUCTURES 


H.  W.  Block 
Z.  Fang 


Abstract  -  Characterization  problems  of  setwise  independence  are  considered  by  introducing  the  set¬ 
wise  dependence  class  SMk(n). 

Two  questions  are  answered:  (i)  Under  which  conditions  are  a  strongly  positive  orthant  dependence 
random  vector  setwise  independent;  (ii)  Under  what  kind  of  dependence  structure  does  a  block  diago¬ 
nal  Covariance  matrix  imply  setwise  independence.  The  techniques  involved  in  dealing  with  the  above 
questions  allow  us  to  answer  some  open  problems  involving  the  SPOD  class  raised  in  Newman  (1984). 
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87-05 


CONCEPTS  OF  SETWISE  DEPENDENCE 


D.  Chhetry 
G.  Kimeldorf 
A.  K.  Sampson 


Abstract  -  This  paper  motivates  and  introduces  a  number  of  new  concepts  of  positive  dependence 
among  sets  of  random  variables.  In  particular,  setwise  positive  upper  (and  lower)  orthant  dependence, 
setwise  association  and  other  related  concepts  are  studied  and  their  relationships  are  explored.  These 
new  concepts  are  applied  to  various  multivariate  normal  distributions,  some  of  which  have  both  posi¬ 
tive  and  negative  covariances.  A  variety  of  new  and  interesting  inequalities  are  also  obtained. 


87-06  A  COST  RELATIONSHIP  BETWEEN  T.  H.  Savits 

AGE  AND  BLOCK  REPLACEMENT  POLICIES 


Abstract  -  The  general  cost  structure  of  a  unit  on  line  is  assumed  to  be  governed  by  a  stochastic  pro¬ 
cess  (R(t),  0  <>  t  <,  Q,  where  R(t)  denotes  the  operating  cost  on  [0,t)  and  £  denotes  the  time  of  an 
unscheduled  (or  unplanned)  replacement  by  a  new  unit  at  a  cost  Cj.  For  an  age  replacement  mainte¬ 
nance  policy,  scheduled  (or  planned)  replacements  occur  whenever  an  operating  unit  reaches  age  T, 
whereas  in  the  block  replacement  case,  scheduled  replacements  occur  every  T  units  of  time.  Such 
scheduled  replacements  cost  C2-  The  expected  long  run  cost  per  unit  time  can  then  be  expressed  in  the 
form  A(T)/E[min(£,  T)]  and  B(T)/T  respectively.  Our  main  result  shows  that 
B(T)  =  j0  t)A(T-x)  dU(x)  where  U  is  the  associated  renewal  function  generated  by  £. 


87-07  NONLINEAR  MODELS  FOR  TIME  SERIES  WITH  ARBITRARY  N.  A.  Langberg 
CONTINUOUS  MARGINAL  PREDICTION  -  THEORY,  D.  S.  Stoffer 

ESTIMATION  AND  PREDICTION 

Abstract  -  In  this  paper  we  present  a  general  method  for  analyzing  stationary  time  series  with  arbi¬ 
trary  continuous  marginal  distributions.  Both  a  nonparametric  approach  (when  the  distribution  is 
unspecified)  and  a  parametric  approach  (when  the  distribution  is  specified)  are  taken.  The  models 
which  we  construct  are  nonlinear  models  that  are  obtained  by  a  monotone  nondecreasing  transforma¬ 
tion  of  a  Gaussian  linear  process.  We  discuss  the  estimation  and  prediction  procedures  for  both  the 
nonparametric  and  parametric  models  and  apply  both  approaches  to  a  real  data  set. 
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87-08  A  MULTIPLE  COMPARISONS  PROCEDURE  FOR  USE  IN  P.  K.  Norwood 

CONJUNCTION  WITH  THE  BENARD-VAN  ELTEREN  TEST  A.  R.  Sampson 

R.  Staum 
K.  McCarroll 


Abstract  -  The  Bemard-van  Elteren  test  statistic  is  a  generalization  of  the  Friedman  test  statistic 
which  permits  multiple  and  unequal  replicates  per  treatment  per  block.  A  multiple  comparisons  pro¬ 
cedure  deveioped  by  Nemenyi,  based  on  the  Friedman  rank  sums,  is  extended  for  the  pairwise  com¬ 
parisons  of  treatments  in  conjunction  with  the  Benard-van  Elteren  test.  This  procedure  is  applied  to  an 
example  arising  from  a  clinical  trial.  A  necessary  and  sufficient  condition  for  use  of  this  multiple  com¬ 
parisons  procedure  is  given,  with  examples  of  designs  for  which  this  condition  is  met.  This  multiple 
comparisons  procedure  is  shown  to  be  identical  to  the  Nemenyi  proceduie  for  the  RCB  design  with 
one  replicate  per  treatment  per  block,  and  to  a  procedure  given  by  Campbell  for  the  BIB  design. 


87-09  PARTIAL  ORDERS  ON  PERMUTATIONS  AND  DEPENDENCE  H.  W.  Block 
ORDERINGS  ON  BIVARIATE  EMPIRICAL  DISTRIBUTIONS  D.  Chhetry 

Z.  Fang 
A.  R.  Sampson 

bt 

Abstract  -  Three  well  known  partial  orderings  and  one  new  one,  designated  'ey  t  =  1,2, 3,4  ate 

defined  on  permutations  of  { l,2....,n}  through  a  unified  approach.  Various  formulations  of  these  par¬ 
tial  orderings  are  also  considered.  With  the  aid  of  these  formulations,  we  show  that  the  four  orderings 
on  permutations  are  equivalent  to  positive  dependence  orderings  defined  over  empirical  distributions 
based  on  rank  data.  In  particular,  we  show  that  the  orderings  bp  b2,  b3  and  b4  are  equivalent,  respec¬ 
tively,  to  more  concordant,  more  row  regression  dependent,  more  column  regression  dependent  and 
more  associated. 


87-10  IMPROVED  RELIABILITY  BOUNDS  FOR  SYSTEMS  Zhaoben  Fang 

WITH  DEPENDENT  COMPONENTS 

Abstract  -  Improved  reliability  bounds  for  systems  with  dependent  components  are  obtained  by  inves¬ 
tigating  the  multivariate  probability  inequalities  under  associated  or  strongly  positive  orthant  depen¬ 
dence  (SPOD)  assumptions.  Using  information  on  correlatedness  for  different  path  sets  (or  cut  sets) 
and  covariances  between  different  components  more  accurate  reliability  bounds  are  obtained.  The  reli¬ 
ability  bounds  for  a  system  can  also  be  improved  by  an  adaptable  probability  inequality  (Theorem  1  in 
paper).  A  numerical  example  is  given  to  show  the  improvements. 
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87-11  METRICS  ON  PERMUTATIONS  BASED  H.  W.  Block 

ON  PARTIAL  ORDERINGS  D.  Chhetry 

Z.  Fang 
A.  R.  Sampson 

Abstract  -  For  two  permutations  i  and  j  of  the  integers  1,  .  .  .  ,  n,  a  number  of  different  metrics 
between  i  and  j  are  studied  in  a  unified  fashion.  Attention  is  focused  on  two  new  metrics  and  on  three 
previously  studied  metrics.  The  new  metrics  are  related  to  the  more  associated  partial  ordering 
between  pairs  of  bivariate  random  vectors.  Results  are  also  given  connecting  these  metrics  to  a 
number  of  previously  considered  partial  orderings  on  the  set  of  all  permutations  of  1,  ...  ,n. 


87-12  A  REVIEW  OF  MULTIVARIATE  INTERPRETATION  OF  A.  R.  Sampson 

CLINICAL  LABORATORY  DATA 


Abstract  -  This  is  a  review  for  the  Journal  of  the  American  Statistical  Associaton  of  Multivari¬ 
ate  Interpretation  of  Clinical  Laboratory  Data  written  by  A.  Albert  and  E.  K.  Harris  (1987). 


88-01  A  DISCOUNTED  COST  RELATIONSHIP  C.  S.  Chen 

T.  H.  Savits 


Abstract  -  In  Savits  (1987)  a  very  general  cost  mechanism  for  a  maintained  system  was  considered. 
There  he  established  a  relationship  between  the  expected  long  run  cost  per  unit  time  for  the  age  and 
block  maintenance  policies.  In  the  present  paper  a  similar  relationship  is  obtained  for  the  expected 
total  a-discounted  cost. 


88-02  EXACT  POWER  AND  SAMPLE  SIZE  CALCULATIONS  C.  Gatsonis 

FOR  OBSERVATIONAL  STUDIES  A.  R.  Sampson 

Abstract  -  In  observational  studies  it  is  often  the  case  that  the  values  of  the  independent  variables 
cannot  be  fixed  in  advance  but  are  themselves  outcomes  of  the  study.  Power  and  sample  size  calcula¬ 
tions  should  take  this  source  of  variability  into  account,  since  they  have  to  be  done  before  the  data 
become  available.  This  article  reviews  the  mathematical  framework  for  studies  in  which  a  multivariate 
normal  distribution  can  be  assumed  jointly  for  the  dependent  and  independent  variables  under  observa¬ 
tion.  Within  this  framework,  exact  power  and  sample  size  calculations  can  be  done  using  a  series 
expansion  for  the  distribution  of  the  sample  multiple  correlation  coefficient  A  new  table  of  exact  sam¬ 
ples  sizes  for  tests  of  level  .05  is  provided  here.  Approximations  to  the  exact  power  are  discussed, 
most  notably  those  of  Cohen,  1977.  A  rigorous  justification  for  Cohen’s  approximations  is  given. 
Comparisons  with  the  exact  answers  show  that  the  approximations  are  quite  accurate  in  many  situa¬ 
tions  of  practical  interest.  More  extensive  tables  and  an  interactive  computer  program  for  both  exact 
power  and  sample  size  computations  can  be  obtained  directly  from  the  authors. 


-  21  - 


88-03  COMPUTATIONAL  ASPECTS  OF  ASSOCIATION  FOR  A.  R.  Sampson 

BIVARIATE  DISCRETE  DISTRIBUTIONS  L.  R.  Whitaker 

Abstract  por  bivariate  discrete  probability  distributions,  P ,  on  the  M  x  N  lattice,  various  aspects 
for  checking  association  (Esary,  Proschan  and  Walkup  ( 1 067))  are  considered.  A  new  algorithm  is 
given  for  verifying  whether  or  not  P  is  associated.  The  efficiency  of  this  algorithm  is  obtained  and 
compared  to  the  efficiency  of  a  simple  algorithm  based  on  the  definition  of  association.  When 
M  sA/  =  4,  for  example,  the  new  algorithm  requires  less  then  3%  of  the  computations  required  for 
the  simple  algorithm.  In  obtaining  these  results  a  new  set  function  Q  is  constructed  from  P ,  on  all 
upper  sets  in  the  lattice.  In  order  to  construct  the  algorithm  to  check  association,  we  define  a  computa¬ 
tionally  important  set  of  extreme  points  and  consider  related  combinatorics. 


88-04  CHARACTERIZING  AND  GENERATING  BIVARIATE  M.  H.  Metry 

EMPIRICAL  RANK  DISTRIBUTIONS  SATISFYING  A.  R.  Sampson 

CERTAIN  POSITIVE  DEPENDENCE  CONCEPTS 

Abstract  -  Let  i  =  (ij,  .  .  .  ,iN)  be  a  permutation  in  SN,  the  set  of  all  permutations  of  {1,  ...  ,N} 
and  define  a  related  bivariate  random  vector  (R,  S)  by  P(R  =  k,  S  =  ik)  =  N-1,  k  =  1,  .  .  .  ,N.  A 
permutation  i  is  said  to  have  a  given  positive  dependence  property  if  the  related  random  vector  has  that 
property.  We  shrw  !  is  a  PQD  permutation  if  i/k  >  (//k)N,  for  all  /  =  1,  .  .  .  ,k  and  all 
k  =  1,  .  .  .  ,N  -  1,  where  ij  t  <  i2,k  <  '  '  '  <^1  k  =  1,  ...  ,N  -  1  are  the  increasing  rearrange¬ 
ments  of  iiM-k+b  iN-k+2’  •  •  •  Jjsj-  Additionally,  we  show  i  =  (1,  .  .  .  ,N)  is  the  only  permutation 
which  is  any  of  TP2.  SI,  RCSI,  LTD  or  RTI.  Explicit  techniques  for  iteratively  characterizing  PQD 
permutations  are  obtained  and  used  to  count  the  number  of  PQD  permutations  in  SN.  Based  on  these 
techniques,  graph  theoretic  methods  are  used  to  generate  all  PQD  permutations  on  S^-  Similar  results 
for  related  dependence  notions  are  also  given. 


88-05  DUAL  QUALITY  CONTROL  FOR  THE  CYTOLYTIC  ASSAY  R.  Day 

J.  Bryant 
T.  Whiteside 

Abstract  -  The  cytolytic  assay  performed  at  the  Pittsburgh  Cancer  Institute  is  described  and  discussed. 
The  desirability  of  a  convenient  method  to  adjust  assay  data  for  systematic  daily  variation  is  docu¬ 
mented,  and  several  such  methods  are  compared.  These  include  naive  adjustments  obtained  by  com¬ 
puting  the  difference  of  the  daily  mean  of  frozen  control  specimens  and  their  corresponding  long-term 
mean,  use  of  an  ad-hoc  shrinking  constant  determined  by  fitting  fresh  data  to  frozen  controls  by  least 
squares,  and  computation  of  adjustments  based  on  empirical  Bayes  methods  and  EM  algorithm. 
Although  final  conclusions  have  yet  to  be  drawn  concerning  the  latter  two  methods,  it  appears  that  the 
naive  estimator  rather  seriously  over-adusts  the  data. 
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88-06  STOCHASTIC  COMPARISONS:  BLOCK  POLICIES  H.  W.  Block 

N.  A.  Langberg 
T.  H.  Savits 

Complete  repair  and  minimal  repair  models  with  a  block  maintenance  policy  are  considered.  Each  of 
these  models  gives  rise  to  a  counting  process  and  these  processes  are  compared  stochastically.  This 
contrasts  with  most  previous  work  on  maintenance  policies  where  only  univariate  marginal  comparis¬ 
ons  were  made.  Also  a  more  general  block  schedule  is  considered  than  is  customary. 


88-07 


OPTIMAL  PRODUCT-TYPE  PROBABILITY  BOUNDS 


H.  YV.  Block 
T.  Costigan 
A.  R.  Sampson 


Abstract  -  The  usual  second  order  product-  type  approximation  to  Pp  =  P(Xi  <  Cj,  X2  ^ 
c2,....,Xn  <  cn)  is  defined  as  P2  =  ITT  P^i  -  1  -  ci  -  i»  Xj  <  Cj)/J^[  P(X;  £  q).  In  this  paper, 

i=2  i=2 

techniques  of  graph  theory  are  used  to  obtain  a  class  of  more  general  product-type  probability  approxi¬ 
mations,  analogous  to  Hunter’s  second  order  Bonferroni-type  bounds.  For  any  spanning  tree,  T,  of  the 
complete  graph  with  vertex  set  V  =  {1,  .  .  .  ,  n},  the  second  order  product-type  approximation 
corresponding  to  T  ,  P2(T)  is  defined  as 

p2(t>  =  n  p(x;  <  c, ,  Xj  <  cp  /  n  p(Xj  s  Ci)j^  - 1 , 

(i.j)eT  i=l 

where  dj  is  the  degree  of  i  in  T. 

It  is  shown  that  Kruskal’s  maximal  spanning  tree  algorithm  applied  to  the  network  function 
f(qj)  =  P(^i  ^  q  ,  Xj  <,  Cj)/[P(Xj  <  q)P  ,  (Xj  <  Cj)]  yields  the  spanning  tree,  T0,  which  maxim¬ 
izes  p2(T)  among  all  spanning  trees. 

For  multivariate  normal  distributions  a  simple  bound  condition  which  guarantees  that  Pn  >  P2(T)  is 
obtained  by  using  Worsley’s  increasing  representation  of  spanning  trees. 

It  is  also  shown  that  when  applicable,  the  second  order  product-type  bound  corresponding  to 
T  ,  P2(T)  ,  is  sharper  than  the  second  order  Bonferroni-type  bound  of  Hunter  and  Worsley.  This 
allows  us  to  obtain  sharper  bounds  than  those  previously  proposed  for  the  detection  of  a  change  point. 
Additional  applications  considered  include  autoregessive  correlation  models,  modeling  communication 
networks,  multiple  comparisons  to  a  control  and  successive  differences  in  population  means. 
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88-08  PRODUCT-TYPE  PROBABILITY  BOUNDS  OF  HIGHER  ORDER  H.  W.  Block 

T.  Costigan 
A.  R.  Sampson 


Abstract  -  Glaz  and  Johnson  introduce  i-th  order  product-type  approximations, 

(3j  ,  i  =  l,  .  .  .  ,n  —  1,  for  Pn  =  P(Xj  <  Ci  ,  X2  ^  c2,  .  .  .  ,X„  <  cn)  and  show  that 

Pn  >  Pn  _  1  ^  Pn  _  2  -  ‘  ‘  -  P2  -  Pi  when  X  is  MTP2.  In  this  article,  it  is  shown  that  Pn  >  P2  = 

n  n-1  n 

np(xl.1sci.1,  x, < C]>i  n  p(x, <  c,)2p,  =  n  P(Xj  <  Cj)  under  weaker  positive  depen- 
i=2  i=2  i=l 

dence  conditions.  For  multivariate  normal  distributions,  these  conditions  reduce  to  COV(Xj  ,  X.)  >  0 
for  1  <  i  <  j  <  n  and  COV(Xj  ,  Xj  |Xj  _  j)  >  0  for  1  <  i  <  j  -  I,  j  =  3,  .  .  .  ,n.  It  is  shown  that 
—  P2  —  Pi  when  group  sequential  analysis  is  performed  for  bivariate  normal  responses. 

Conditions  for  Pn  >  P3  >  P2  -  Pi  are  a'so  derived. 

Bound  conditions  are  also  obtained  which  insure  that  product-type  approximations  are  nested  lower 
bounds  to  upper  orthant  probabilities  P(X[  >  Cj,  .  .  .  ,Xn  >  cn).  It  is  shown  that  these  conditions  are 
satisfied  for  the  multivariate  exponential  distribution  of  Marshall  and  Olkin. 


88-09  POSITIVE  DEPENDENCE  CONCEPTS  FOR  MULTIVARIATE  M.  H.  Nletry 

EMPIRICAL  RANK  DISTRIBUTIONS  A.  R.  Sampson 


Abstract  -  The  study  of  positive  dependence  notions  for  bivariate  empirical  rank  distributions  can  be 
viewed  mathematically  as  the  study  of  partial  orderings  on  the  set,  Sjsj,  of  permutations  of 
{1,  .  .  .  ,N}.  This  paper  extends  many  of  the  previously  obtained  bivariate  results  to  cases  of  mul¬ 
tivariate  empirical  rank  distributions,  with  primary  focus  on  the  trivariate  case.  Mathematically,  this 
can  be  reduced  to  studying  notions  of  positive  dependence  on  the  space  SN  x  SN  via  the  concept  of 
related  random  vectors,  where  a  tnvariate  random  variable  (R  ,  S  ,  T)  is  said  to  be  related  to 
(i  ,  j)  e  Sfj  x  Sn  if  P(R  =  k  ,  S  =  ik  ,  T  =  jk)  =  N-1,  for  all  k  =  1,  .  .  .  ,N.  Characterizations 
of  the  positive  dependence  properties:  PUOD,  PLOD  and  POD  are  obtained  in  terms  of  a  newly 
defined  notion  of  a  relative  increasing  rearrangement  Also,  characterization  of  the  positive  dependence 
orderings:  more  PUOD,  more  PLOD  and  more  POD  are  obtained  by  comparing  certain  relative  rear¬ 
rangements.  By  our  results  and  certain  bivariate  results,  a  new  partial  ordering  describing  concordance 
on  Sjq  x  Sn  is  introduced.  Its  connection  with  the  orderings  more  PUOD,  more  PLOD,  and  more 
POD  is  established. 
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88-10  SOME  NEW  PARTIAL  ORDERINGS  ON  SD  and  H.W.  Block 

Sn  x  SD  D.  Chhetry 

Z.  Fang 
A.R.  Sampson 

Abstract  -  A  new  class  of  partial  orderings  on  Sn,  the  set  of  permutations  of  { 1,  ....  n},  is  given. 
Each  of  these  partial  orderings  is  shown  to  be  a  subordenng  of  a  recently  described  partial  ordering  on 
Sn  (Block,  Chhetry  and  Sampson  (1987))  which  is  related  to  Schriever’s  (Annlas  of  Staistics,  1987) 
more  associated  ordering  on  bivariate  distributions.  Also  given  is  an  extension  of  three  known  partial 
orders  on  Sn  to  partial  orders  on  Sn  x  Sn.  These  extensions  facilitate  the  study  of  functions  from 
S„  x  S0  into  R,  which  preserve  these  partial  orderings,  thereby,  providing  a  methodology  for  extend¬ 
ing  the  notion  of  arrangement  increasing  functions. 


89-01  LIMITING  DISTRIBUTIONS  FOR  STATISTICS  BASED  S.K.  Bar- Lev 

ON  DIFFERENT  SAMPLES  N.A.  Langberg 

B.  Reiser 


Abstract  -  Let  and  be  two  k-dimensional  statistics  based  on  independent  samples  of  sizes  m 
and  n,  respectively,  from  two  populations,  and  let  f  be  a  vector  valued  function  defined  on  the  2k-th 
Euclidean  space.  This  paper  deals  with  the  limiting  distribution  function  of  f(Am,Bn),  suitably  normal¬ 
ized,  as  m  and  n  tend  to  infinity  separately.  Applications  to  estimation  and  hypothesis  testing  problems 
are  given.  The  connection  between  these  results  and  those  previously  obtained  in  the  literature  is  dis¬ 
cussed. 


89-02  SEQUENCING  OF  EXPERIMENTS  FOR  LINEAR  AND  J.  Bryant 

QUADRATIC  TIME  EFFECTS  G.  Constantine 


Abstract  -  Run  orders  that  are  optimal  against  a  linear  time  trend  and  efficient  against  quadratic  trend 
are  constructed.  The  emphasis  is  on  settings  that  allow  d  replications  on  each  of  m  treatments.  A 
description  of  the  linearly  optimal  and  quadratically  most  efficient  sequences  is  also  given.  Explicit 
formulae  are  found,  along  with  asymptotic  bounds,  for  the  efficiency  of  the  sequences  constructed. 


89-03  OPTiMAL  REPLACEMENT  FOR  A  GENERAL  MAINTENANCE  MODEL  C.S.  Chen 

T.H.  Savits 


Abstract  -  We  continue  the  study  of  our  general  cost  structure  for  a  maintenance  system.  Here  we 
focus  on  the  optimizanon  questions.  The  notion  of  a  marginal  cost  function  is  rigorously  formulated 
and  its  utility  investigated.  Various  applications  are  considered  including  a  new  model  in  which 
minimal  repairs  are  performed  as  long  as  the  total  repair  costs  do  not  exceed  a  fixed  amount. 
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89-04  A  NOTE  ON  CHARACTERIZATION  OF  A  MULTIVARIATE  T.T.  Nguyen 

STABLE  DISTRIBUTION  A.R.  Sampson 

Abstract  -  Several  characterization  of  multivariate  stable  distributions  together  with  a  characterization 
of  multivariate  normal  distributions  and  multivariate  stable  distributions  with  Cauchy  marginals  are 
given.  These  are  related  to  some  standard  characterizations  of  Marcinkiewicz. 


89-05  SECOND  ORDER  BONFERRONI-TYPE,  PRODUCT-TYPE  AND  H.W.  Block 
AND  SETWISE  PROBABILITY  INEQUALITIES  T.  Costigan 

A.R.  Sampson 

Abstract:  This  review  paper  considers  a  variety  of  lower  bounds  to  Prob[X[  <  c,  .  .  .  ,  XD  <  c] 
involving  one  and  two  dimensional  marginal  distributions.  Some  of  these  bounds,  e.g.,  Bonferrom- 
type,  do  not  require  conditions  on  the  random  variables.  Others  of  these  inequalities,  e.g.,  product- 
type,  require  positive  dependence  conditions  on  Xj,  .  .  .  ,  Xn  for  the  inequalities  to  hold.  Because  all 
of  the  two  dimensional  bounds  depend  on  the  labeling  order  of  the  random  variables,  various 
permutation-optimized  versions  of  the  bounds  are  described.  Relationships  among  the  various  bounds 
are  also  considered. 


89-06  REPAIR  POLICIES  AND  STOCHASTIC  ORDER  T.H.  Savits 


Abstract  -  This  paper  reviews  the  role  of  stochastic  order  as  it  relates  to  the  study  of  maintained  sys¬ 
tems  in  reliability  theory.  The  classical  univariate  comparison  results  for  the  age  and  block  replace¬ 
ment  policies  are  presented.  Extensions  to  stochastic  comparison  of  processes  and  recent  generaliza¬ 
tions  of  age  and  block  policies  are  also  discussed. 


89-07 


REPAIR  REPLACEMENT  POLICIES 


H.W.  Block 
N.A.  Langberg 
T.H.  Savits 


Abstract  -  In  this  paper  we  introduce  the  concept  of  repair  replacement.  Repair  replacement  is  a 
maintenance  policy  in  which  items  are  preventively  maintained  when  a  certain  time  has  elapsed  since 
their  last  repair.  This  differs  from  age  replacement  where  a  certain  amount  of  time  has  elapsed  since 
the  last  replacement.  If  the  last  repair  was  a  complete  repair,  repair  replacement  is  essentially  the  same 
as  age  replacement.  It  is  in  the  case  of  minimal  repair  that  these  two  policies  differ.  We  then  make 
comparison  between  various  types  of  policies  in  order  to  determine  when  and  under  which  condiuon 
one  type  of  policy  is  better  than  another. 
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